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1 Introduction 



String theory is the best candidate for the unified theory of interactions. Dirich- 
let branes Q] on which fundamental open string boundaries are attached are key 
objects to understand non-perturbative effects of string theory. That is, in the 
type I1B string theory, the D-string is a strong-weak coupling dual object of 
the fundamental string. Furthermore a matrix model for 11B string theory was 
presented by IKKT §]. They propose that the model is a non-perturbative def- 
inition of string theory. However, clear correspondences between string theory 
and its matrix model are not found yet. Of course in Ref some correspon- 
dences between them were presented. Firstly, Green Schwarz action in Schild 
gauge can be related to the action of the matrix model in the case of the infi- 
nite matrix size. Next, the potential between two static D-strings in the matrix 
model was calculated. In Ref jj| the authors calculated the scattering ampli- 
tude of two D-strings based on the matrix model and they compared it with the 
result j| in string theory. They set the non-zero values of all commutators of 
the coordinate matrices to a 1 x constant. However in Ref H, the authors insist 
that the value of the commutators in the matrix model F is the inverse of the 
field strength F on the D-brane world volume in string theory in the large F 
approximation. In this context, the authors in Ref Q compared the potentials 
between two static D-branes in string theory and the IIB matrix model. 
We consider Dp-branes with non-zero field strength in IIB string theory and in 
its matrix model (IKKT model). We calculate the amplitude of the scattering 
of two Dp-branes and compare it with the one Q of the IIB matrix model. 
Since the relative velocity (2v) is non zero, the configuration is not BPS. When 
v — > 0, the amplitude is expected to vanish. We find that the amplitudes in the 
two models are equal in three independent limits if we identify the respective 
field strengths suitably. Namely, we find how the "D-brane" configurations in 
the matrix model are translated to those in string theory. 

In section 2, we calculate the amplitude of two Dp-branes with the field 
strength on their world volumes in bosonic sector. First, we explain the con- 
figuration of two D-branes in the flat space-time. Two Dp-branes are extended 
infinitely in the world volume direction. They have a constant relative velocity 
to each other and flat world volume. They are located in spatially parallel to 
each other. In this situation we consider the fundamental open string whose 
ends are attached on the two D-branes. Then we calculate the one loop ampli- 
tude (i.e. cylinder amplitude) of open string. When the equations of motion 
are chosen to be free field type, the boundary condition of open string is com- 
plex because of the non-zero field strength F and the non-zero relative velocity 
2v. Previously, v — case was calculated in ||, 0, F — case was also done 
in || and related case in IIA string theory was presented in ||] . We formulate a 
general mode expansion of the string coordinates in the flat world volume case, 
quantize them canonically, and construct the Virasoro algebra. 

In section 3, we consider superstring based on the NSR formalism. In this 
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paper, the boundary condition for the world sheet fermion is defined by the 
super transformation of the bosonic one. Then we perform a calculation using 
the proper time representation of the amplitude. To compare the result with 
that of the matrix model, we choose the leading term in three independent ex- 
pansions in terms of parameters, in the long distance (b ^> a') case where b 
is the distance between two D-branes, in the low velocity case, or in the large 
gauge field strength case. 

In section 4, we consider the amplitude in the IKKT matrix model with a 
brief review in our notations. 

In section 5, we compare the amplitudes in the two models and identify the 
field strength F in string theory with the value F M of the commutator in the 
matrix model, in such a way that det (77 + F)F M = — 1 where 77 is the target 
space metric. Then we find that the amplitudes in the two models are precisely 
the same when b is large for arbitrary v, F, v is small for arbitrary 6, F, or F is 
large for arbitrary b, v. when the distance between the two Dp-branes is large 
for arbitrary velocity v. 

In section 6, there are conclusions and discussions. In Appendices, we ex- 
plain a general method to study the open string boundary condition which is 
any linear combinations of Neumann and Dirichlet boundary conditions. Some 
detailed calculations are also explained there. 

1.1 physical situation 




First, we consider two flat identical D-branes with the same world volume di- 
mensions in the D = 10 dimensional Minkowski space-time whose coordinates 
are written as X°, • • • , X D ~ X . We restrict the world volume dimension: (p + 1) 
to be even because we consider the IIB string theory. We locate them spatially 
parallel to each other in the D dimensional Minkowski space-time. We call the 
common directions of the two Dp-brancs 1, • • -p, where the world volumes of 
the Dp-branes are infinitely extended. We set two Dp-branes to have relative 
velocity, 2v, in the (p + l)-th direction. We set the coordinates' value of the 
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p + 2, ■ ■■ ,(D — 2)-th directions of the two Dp-branes to be and locate them 
separately in the (D — l)-th direction by the distance b. We restrict v and b to 
be positive constant in this paper. 



1.2 correspondence between string theory and matrix model 



Tree 



/Classical 



D-b 



\backg 



rane 
round/ 




Matrix 
Picture 



World Sheet 
Picture 



2-Body Planar graph Interacting branes 



Assuming that two D-branes interact by exchanging closed fundamental strings 
and any effect of two D-branes can be treated as open string boundary condi- 
tions, we calculate the amplitude of the fundamental string diagram with two 
boundaries, especially the cylinder amplitude. Since what is calculated in the 
matrix model is one loop diagram which is planar, we expect that the cylinder 
amplitude agrees with the sum of planar diagrams in the IIB matrix model. The 
cylinder amplitude can be thought of the one loop amplitude of an open string 
whose boundaries are fixed in the D-branes, which we are considering now. 
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2 D]?-brane-Dj)-brane scattering for bosonic part 



We would like to calculate the scattering amplitude, that is the effective action. 
To calculate the amplitude, we have to obtain the propagators in space-time, 
in another word, the normal ordered Hamiltonian with the correct zero point 
energy. Therefore, in the following part in this section we define a normal 
ordering, construct Virasoro operators, and find a physical state condition. Let 
us start with the following world sheet action of open string. 



Si 



1 



bosonic 



Ana 

1 



2ira' 



dcrdTd a X^d a X tJ 



dTAld T X" 



2ira' 



(2.1) 



where and A™ are the gauge potentials on the D-brane world volume to 
which the a — and a — ir end of open string are attached, respectively, lj 
We can rewrite the action as follows. 



Sbosonic = {terms whose deformation vanishes with free equations of motion} 

~l~ ^boundary 5 



where 



5', 



boundary 



47TQ 



■jX^ (d a X^ - 2d T Al) 



(2.2) 



To get free equations of motion, we require 



5S-. 



boundary 



= o 



'In this paper, the indices run as follows: 



fi, v, p, a = 
a,/3,7 = 0, ■ 
A,<5 = 2, 
a, 6 = 0,1 
i,j = 2,-- 
k = p + 2, 
1,1' = I,-- 
m, n 6 Z 
r, s £ Z or 



■ D — I 



D-l 



■ ■ D — 

' 2 



(2.3) 



in the D=10 dimensional space-time 
in the p-brane world volume 

in the pure Neumann and pure Dirichlet directions 

the meaning depends on the situations 

in the mode directions 

in the pure Dirichlet direction 

for the field strength 

for bosonic modes and b, c ghosts 

for fermionic modes and /3, 7 ghosts 
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Namely, 



6X» (d a X,-2d T Al)\ a=0 
5X» - 2d T Al 



V>/ I (7—7T 



. 



(2.4) 



2.1 p — 1 case 



In this subsection, we consider p = 1 case to understand how to treat the non 
trivial boundary conditions. 

First, imagine a pure (with no gauge field) D-string moving in the 2-nd direction 
with the constant velocity v. Its coordinates are restricted as 



it is equivalent to 



X 2 = vX° 



d T (x 2 - vx°) = o . 



(2.5) 



(2.6) 



This is a Dirichlct boundary condition. 

Second, in order to introduce gauge field on^the world sheet of the D-string, 
look at it from the proper time coordinate: X in the following way. 



so that 



where 



d a x» — d a x» = h» v d a x v , 



d T X 2 = 



A H _ 
1 '-v v — 











1 

1 







V 







(2.7) 
(2.8) 

(2.9) 



1 / 



The boundary condition in this coordinate is, because of Lorentz invariance, 

6X" (drXp - 2d T A *) = , (2.10) 
on each end. To consider D-string, we impose the following conditions. 

5X°,SX 1 ^Q,SX 2 = --- = 6X 9 = 0. (2.11) 

Therefore 



d a X - 2d T A%* = 
d a X 1 - 2d T A° 



0,TT 







(2.12) 



d T x< 



d T X 9 = 



G 



We consider the following constant field strength of the gauge field. 



a/3 



o A 
-fi o 







(2.13) 



where .F Q| g is the same in each brane since we are considering identical branes. 
Now, the boundary condition is 



8 a X + fidrX 1 = 

d T X 2 = ■■■= d T X 9 = 

on each end. 

We can rewrite the boundary condition as 



(2.14) 




By using the fact that 



X 1 — X\, X a — —Xq 



because of Minkowski space-time, we can rewrite eq. (2.15) again as I 



dX L = M dX R , M := 



/ i+/i 2 

/ 1 71 



2/l 

i-/r t=7? 

2/i 

T=7? i=7? 
o 



. (2.15) 



(2.16) 



(2.17) 



where we decomposed the coordinate functions by using the equations of motion 

X(r,a) =X L (a+) + X R (a-) , a ± := r ± a . (2.18) 

Then finally, we can get the following formula for the boundary condition in the 
original coordinate flame: 



dX L = M v dX R ■ ■ ■ a = 
dX L = M_ v dX R ■ ■ ■ a = tt 



(2.19) 



4 We sometimes omit indices. Eq. (2.17) means 

CT — 0,7T 



<T — 0,7T 
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2.2 mode expansions for general boundary conditions 

We can treat general boundary condition with two Dp-branes as follows. 



( 3X l {t) = MdX_ R { T ) •••o- = 

\ dX L (r + tt) = MdX R (r - tt) • ■ • ct = tt 



(2.20) 



where M and M are any SO(D — 1, 1) matrices. In this section, we proceed in 
this general case for a while, and sometimes come back to our special case: two 
identical parallel Dp-branes' scattering. Next, we have to perform mode expan- 
sions of the coordinate functions which satisfy the above boundary condition. 
The result is as follows. 



X(t,<t) = x + 



E 









M [ - 




dpe 


Jo 


Jo 





-i(m-\-iE) p 



where 



E := — InM^M 
2tt 



(2.21) 



(2.22) 



is restricted to be diagonalizable. x is the integration constant for r and cr, 
which may depend on a m 's. 



2.3 quantization 

We perform the canonical quantization by applying the canonical commutation 
relations to the coordinates and their conjugate momenta. The only non-zero 
commutation relations are, 



[d T X"(T, a), I"(t, a')} = -2i:ia'S(a - J)rf ', (0 < a,a' < n) , (2.23) 



where 



X] = diag(-l, 1, • • • , 1) 



We can get the commutation relations between modes from eq. ( 2.23] ) . 
The results are B □ 



\X. (Xr. 



(m + iE) r)5 m+n , 



ma 



(M - M~ x ) T) 



(2.24) 

(2.25) 
(2.26) 

(2.27) 



5 See Appendix A. 

6 For any Lorentz vectors and b" , matrix [a, b] is denned by its elements to be [a^ ,&"]. 
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Note that the component of the matrix M in the pure Neumann or pure Dirichlct 



direction are 1 
convenient form 



respectively. We can rewrite these relations as a more 



m + iE\i 
m 





N 22 

Nn coth 7tEh7]t n-^ii 
+JViH 7 T 11 (ST)* 1 -(ST)ii 7) T 11 JV 1 t 1 




where 



:= T a r < 



Sx 



2 r 

gZ 



C i) 



T~ lr qT~ 



(2.28) 
(2.29) 
(2.30) 

(2.31) 
(2.32) 

(2.33) 



In equations ( 2.28| )-( |2.33| ), matrices T, S, N and C m are defined by M and E 
in the appendix B. In the case of two identical parallel Dp-brane scattering, the 
commutation relations are as follows 



On 




/ 2o'diagf-, 1 + m 2 , •••,! + m* P+ i 



9 ■= 



l-ft(l-v»), i,j = 2,---,p + l, 



where v is the velocity of the D-branes in the (p+l)-th direction and 1112, 
are the 2x2 matrices in the following form. 



mi 



l + /i 2 2fj 

2/1 1+/? 



TO, 




; = 2,---, 



p+1 



(2.34) 

• , mp+2 
2 

(2.35) 



7 See Appendix B. We would like to have complete pairs of canonical variables and their 
conjugate momenta with one to one correspondence, by transforming (x, a m ) — + (x, ct m ). 
8 In this section and appendix B, we write any matrices A as: 



/An Ai 
V Mi A 2 



2 

A22 J 



where the dimension of the square matrix An is defined to be equal to the rank of E. 
®t means transpose. 
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In eq. (2.35) /i, • • • , f p+t are the components of the U(l) gauge field strength: 

0= P 



0=0 

/ f x 

-h o 



F a (3 '■ = 



\ 

fp+i 
-fp+l o j 



a/3 



defined in the proper flame of branes, as in the subsection 2.1 



(2.36) 



2.4 definition of the vacuum 

We consider two identical parallel Dp-branes scattering here. To define normal 
ordered operators, for instance, Virasoro operators, we identify the creation 
and the annihilation operators with respect to the vacuum. Since e ^ 0, the 
commutation relations for the modes include the following unusual relation as 
compared with the case of pure Neumann or pure Dirichlet. 



Therefore we classify the modes as 



a m>0 ' 
a m<0 ' 



a.Q ■ ■ ■ annihilation operators 
• • • creation operators 



That is, the vacuum 10 > is defined as 



al\0 >=< 0\a° = a£ >0 |0 >=< 0|< <0 = 



(2.37) 



(2.38) 



(2.39) 



Since exchange the role of &q and a\ correspond to exchange the sign of the 
imaginary part of the Hamiltonian, we choose the role as above in order to the 
imaginary part to be negative. (Recall that v is assumed to be positive.) 

2.5 Virasoro algebra 

We define the energy momentum tensor as 

/ 1 ^^bosonic 



T ab := -4ira 



Vh Sh ab 



(2.40) 



where Sbosomc is the action in eq. (2.1) before the conformal gauge is chosen. 
We introduce the holomorphic energy momentum tensor, 



T{z) 



1 



2a'z 2 



(T aa + T 01 ) 



(2.41) 
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where 



If we define 



then 



z := e la± . (2.42) 



T(z) := T(z) - ^i±^ , (2.43) 



T(z)T(w) = , i + 2T{W \ 9 + + regular (2.44) 

(z — wp (z — w) z — w 

This is the same form with the simpler boundary conditions. 
Virasoro operators are defined by 

L m := / ^z m+1 f{z) ■ (2.45) 

J\z\=l 2m 

Now we can find the Virasoro algebra: 



Lm-, L r 



f ^ w n + l Res z m+1 T(z)T(w) 
J 2ni z=w 



= (m - n)L m+n + ^m(m 2 - l)S m+n . (2.46) 



Note that the physical state condition by Lq is different from the usual one, 
namely 



(l - cj |phys) 

t ie(l + ie) , 
: a_ n r]a n : c |phys) 




= , (2.47) 

where : : means the normal ordering defined in the previous subsection and c is 
the intercept, which is 1 in the bosonic open string theory. Now we can calculate 
the one loop scattering amplitude for bosonic string. The formula is Jh} 



^bosonic = - In dot ^ (l -cj = i Tr In [h - cj 

= _I r^Tre-t^o-c) (2 48) 

2 Jo t 

However we do not calculate this and we consider the superstring case in the 
next section. 



3 See Appendix C. 
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3 D]?-brane-D]?-brane scattering in superstring 
theory 

In this section we consider supersymmetric case. Let us introduce super partners 
for dynamical variables X in the previous section. We choose RNS formalism 
in this paper and introduce fermion variable ip in the world sheet. In the super 
conformal gauge, the total action is til 



SstTl 



1 



Aira 
1 



7 / dodr [d a x»d a x vi - i^p a d a ^] 



. drF" 
Ana 1 I 

(ghost) 



a=0 



(3.1) 



where 



rpv 



F 




pa 



where ± correspond to a — and a = 7T, respectively and 

p 



A p - 



fj,=0 




H=p+1 


1 




1 




VI— v 2 




Vl-v 2 




— V 






1 


VI— u 2 


VI— v 2 





1 J 



(3.2) 



(3.3) 



In eq. (3.2), F defined in eq. (2.36) is a field strength on D-branes world volume 
and is set to be constant. Since we are considering identical branes, F is the 
same in each brane. A„ is a Lorentz boost matrix with the velocity v in the 
(p + l)-th direction. 
Since the spinor is written as follows 



4>R 
IpL 



the world sheet global supersymmetry is expressed as 

j 6 susy X» = ie susy (tpl + ^ a R f 

\ 5 susy i>7 = -2e susy e^~ a ^d ± X» 



(3.4) 



(3.5) 



Our convention of two dimensional Dirac matrix is 



<M° o ).->-(! ;) 
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where e SU sy is independent of t,ct and 

VA:=e^ CT ±^ •••( a = ° NSsector . (3.6) 
r r L a = 1 R sector 

Due to the existence of the boundary (i.e. D-brane), there is only one parameter 
for supersymmetry 

3.1 boundary condition 

As in bosonic case, in order to get free equations of motion, we require 

SX {r)d l X - F"d X) + (tjpV + F°pV) = . ( 3 - 7 ) 
or represent equivalently as 

SX[( v + F v )dX L -( v -F v )dX R ] = 0, (3.8) 

^(77 + ^)^-^(77-^)^ = 0. (3.9) 

In the proper time coordinate frame of D-branes, we impose-the following bound- 
ary conditions for bosons and fermions at the both ends. tHl 



• Neumann direction 

• Dirichlct direction 

(a = 0) 
(a = tt) 



(3.10) 



(3.11) 



(rj + F)dX L = {r,-F)dX R 

(77+ F)r L =^(v- F)r R 

dx L = -dx R 

1 n = -n 

In the original frame, that is 

8X L = M v dX R 
?Al = M v ip R 
dX L = M_ V 8X R 

i>L = (-) 1 - a M^ I 

12 We determine boundary conditions for fermions as follows. If the boundary condition for 
bosons is 

(d+ -Md-)X = , 

from cq. ( fc.jjj ) as in the section 2, then we transform this by the supersymmetry and get the 
following equation. 

(d+-Md.)(r L + r R ) = o. 

By using the equations of motion for fermions: d±ip a R = 0, we get d T (i>l - Mip%) = 0. 

L 

Then we impose the following boundary condition for fermions. 

V>£ = M^ a R . 

Eq. (13 . 91) determines the relation between Sip^ and Sipn- 
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In eq. (3.11), 



M v := A^MqA. 



(3.12) 



mi 



Mn := 



rrip+i 



(3.13) 



Vo 2 -i/ 

where mi, • • • , mp+i are 2x2 matrices which are defined in the previous section. 

2 

3.2 mode expansion 

We can find the mode expansions for X and ip to satisfy the above boundary 
condition as follows. 



X(t,ct) = x + \/ T ^ 









M v f -+ 


f 


dpe 


Jo 


Jo 





a'M v J2 ^ %{r+lE)a+ dr , 



where 



E:= -L In M^M-v 



(3.15) 
(3.16) 

(3.17) 



which is diagonalizable. 
3.3 quantization 

We can find that the commutation relations for fermionic modes d are the usual 
type: 



{d r ,d s } = r]5 r+s , 
and bosonic commutators is as in the previous section: 
[a m ,a n ] = (m + iE) i]5 m+n 
[x,a m ] = (1 + M v )r) 

[x,x] = (M v - m- 1 ) v 



(3.18) 



(3.19) 



Eq. ( 3. IS ) is rewritten as eq. ( 2.34| ). The ghost part is the same as usual 
because the boundary condition does not change it. 
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3.4 Virasoro algebra 

The energy momentum tensor is 

1 



T(z) 



d-xi*d-X u 



OS Oil 

- E 

2 ^ 

■r.s^'Zj — 

(ghost) . 



crmion 



(z) + Tghost(^) 



(ghost) 



d^r)(s + iE)d s 



If we define a shifted energy momentum tensor as 

ie(l + ie) aD 



T(z) 



T{z) 



„, „ 1 / a_D ie 



2z 2 
i 

16 



16z 2 



then we get 







T(z)7» 



2T(w) d w T{w) 



(z — w) 4 (z — w) 2 z — w 
(regular + ghost contribution) 



This form is of usual type. 



3.4.1 physical state condition 

Virasoro operators are defined by 



— z m+1 T{z) . 
|*|=i 2m 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



Note that the physical state condition by Lq is different from the usual one 
because energy momentum tensor T which satisfies the correct operator product 
expansions is shifted from T as eq. (3.21). That is, 

1 - a s 



Lo--- 



|phys) = . 



(3.24) 



'See appendix C. 
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We further define 



L - L *-2- — 



L« + J Ld + L bc + ^ 7 _| + l_^ ! (3 25) 



Lq := - ^2 '■ a -nmn ■ , (3.26) 

oo 

L o := £ rf -^( r + + 1 4*7^0 , (3.27) 

OO 

L£ c := ^ n (&_„c„ + c_„&„) , (3.28) 

71=1 

OO 

£q 7 : = £ r(/?_ r7r - 7 _ r /? r ) , (3.29) 
where 6, c, /3, and 7 are the standard ghosts. The normal ordering procedure 



in eq. (3.26), is found by eq. (|2.38j). 



3.5 scattering amplitude 

We calculate the scattering amplitude for two Dp-branes, namely, the one loop 
amplitude of open string whose ends are attached on the Dp-brancs. The one 
loop amplitude can be expressed by the following proper time (t) integral of the 
exponentiated Hamiltonian, which annihilate the physical state. It is gotten in 
the previous subsection, that is Lq with the correct constant shift. The shift can 
be found by imposing the physical state condition in the previous subsection. 
Therefore we use Lq as Hamiltonian. The integral region is from to 00 because 
of cylinder Since we are choosing the NSR formalism, we have to sum over 
the spin structures which have 2x2 = 4 combinations because of the cylinder 
amplitude. The formula of the amplitude is |10| 

String = - In det - 3 L a = - Tr In L a = - - / - Tr e^ tL « 

A 2 In t 



= "jf T £ C(l)^-re---, (3.30) 

J{ > a,6=0,l 

where F — F d + F@ n is the fermion number operator, the summation is done 
over the spin structures which consist of GSO projection (b = 0, 1) and the 
summation of NS (a = 0) and R (a = 1) sectors. C(f) are their weights: 

°Q-°(?)- c Q-i- 
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" Tr " means 

NSR 



Tr = 2 x Tr x Tr x Tr x Tr , (3.32) 

NSR a d b,c 0,-y 

where the factor 2 is caused by the sum of the configurations of the exchanged 
ends of open string. Therefore 

Tr (-) bF q L « = 2 x q-i^+'-^r) 



NSR 



x Trq L o x Trq L °' C 

a b,c 

x Tr(-) 6F V° x Tr(-) bFf) '\ L o~' , (3.33) 

d /3,7 

where 

q := e _7r * . (3.34) 
We calculate these traces as follows. 



Trg ° = g 2 ^fc= P +2 xTrg 5 ^.=2 a ° x Tr q 1 frequency 

/"CO P+ 1 



term 



oo 

(1 - Y[ (1 - <f )~ (I? ~ 2) |1 - <f- ie |" 2 , (3.35) 



TrgV* = ^(l-g") 2 , (3.36) 

b.c 

n=l 

Tr(-) bFd q L ?> = (l_ a 6) (Vcos^V 



X i_ s \(fl-2) 



X 

n=l 



JJ (l + (-)V" V )' '' l + {-) h q n -^- u , (3.37) 

a=l 

oo _„ 

Tr(-)^V°^ = Jj(i + (_)V-^)- . (3.38) 

n—1 

In the above equations, the zero mode integral can be done by using the 
relation between p 1 and a 1 obtained in the previous section. .The result can be 
expressed as follows by using Jacobi's theta functions: 9 (?) £3, 



Tr^o+^T) = V p gf[(l + f?)(4n 2 a't) 



-I —£-t 



a, b.c 

1=1 



Vp is uic space volume 01 a jj/j-urane. oee appencux o anu sec also 
15 The definition and some properties are in an appendix in Ref. Jll| . 
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(i)( 



ei I it_ 
2 I 2 



Tr <-) bF 'a( L o+ L o-~') 
(1,(9,7 



(3.39) 
(3.40) 



Therefore we can find the scattering amplitude cq. ( p.30[) as 

p+i 



^string 



i(27r) 3 ^n( 1 + ^) 



;=2 



f 



(4-7T a't) 2 e 4m 



E c 



::)(o|f) 3 ^c)(f if) 



..fei V67 9(1/ (o|f) fl(i)(f If) 

Using Jacobi identity: 

E^:)K:)(o|f) 3 c)(fif)^(;)(!if) 4 . 

a, 6=0,1 

and the property of the theta function under the modular transformation: 



(3.41) 



(3.42) 



^O(flf) = v^^^H-J^It) 

= (f)^(l)'(o|f) 



(3.43) 



we can further simplify the amplitude eq. (3.41) as 



v P gi[(i+m- j(^*> 

1=2 ^° 



. _ P b 

t) 2 e 4, 



(TTi) 3 6 


G)(f 


f) 4 


*G)>lf) 3( 




If) 



(3.44) 
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4 D]?-brane-D]?-brane scattering in IIB matrix 
model 



In this section, with our notation, we briefly review the calculation of the brane 
scattering amplitude in the IIB matrix model |J. In this model, dynamical 
variables are "coordinates" represented by N x N matrices and N is assumed 
to be infinitely large. 
The action is 



a Tr 



±[x lt ,x v ] 2 -±$r'>[x ll ,ti>] ) - -IX 



The classical equations of motions for X M from this action are 

[X> t ,[X li ,X v ]] = 0, 



(4.1) 



(4.2) 



which has the solutions that may be representing D-branes. In our situation: 
two identical parallel p-branes moving with relative velocity 2v, the following 
classical configuration is suitable H ■ 



x (l) 





X (2) 



^ = 0, 



(4.3) 



where Xm are x -y matrices and their indices I = 1, 2 mean the first and 
second branes. 

First, in the proper time frame (written as X) for each brane, we set 



X (i), X , 



= 2ma'F 



I pM a/3-. 



F 



A I 



( o si 
1 -fi 1 o 



V 



. n , all other commutator = , 
\ 





_ fM 
J p+1 



J p+1 

2 

J 



detF M ^ 



(4.4) 



(4.5) 



in order to describe p-branes. 

Since we are consideri ng identical branes, F M is independent of / = 1,2. The 
right hand side of eq. (4.4) is proportional to the identical matrix. To satisfy the 
cq. (4.4), we use 2i— pairs of matrices (qi , pi) with large enough dimensions. 



' -<>22-2 
X (I) 



X 



21-1 
(I) 

xr +1 



D-2 



X 



D-l 



X 



[luPv] = iSw , 1,1' = 1, 



P + 1 



(4.6) 



(i) 



±1 
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Next, for the first and second branes, we perform the boost with velocity v and 
— v along the (p + l)-th axis, respectively as 



±v \ / j^O 



(I) _ VT^ VT^ (I) 
±v 1 vP +1 



where ± correspond to / = 1, 2. We can get the following solution: 



(4.7) 



\ ^VT^ 



qi o \ _ p+i 



JSf2i-l _ / ^ ^ 7 = 1... £±1 

-\l**aj l I o pj ' X ' ' 2 (48) 



go 



— u 



^ p+2 ...D-2 _ q 

v-D-l _ / 3 
-§ 



where 2i> is the relative velocity between the branes. 

The result of the calculation for the effective action (scattering amplitude) in 
the one loop level is [[| 



2 />00 l, 

A matlix = -W p t\fr ■ / T( 47rV *)" S 

*=2 70 4 

. 4 //f sinhe A/ A 
1 Sm ( ' 2 *j 



e inc.' 1 



ff* cosh e M sin sinh e M t) ' 



(4.9) 



where 



tanhe M :=v, (4.10) 
and the fact that the matrix size N is large enough is used. 
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5 Comparing string theory and matrix model 

We compare the two scattering amplitudes in the IIB string theory and the IIB 
matrix model, respectively. They can be rewritten as follows: 



^string 



2 Vp 



n (>? + !) 



n 



1 - ft 1 - V 

. . 4 

sinh 4^ 



dt . t o . . _ P _ b 2 + 

, — (An a t) 2 e w 1 
o * 



n=l 1 



/ sinh 7Tg \ 



2 ' 



2lVl - f2(l - V 2 ) , . 



1 matrix 



2 p p+1 



n 



Z=2 



J t 



dt 



— (47r a't) 2 e i^ 7 



, sinhe M = 



=i 1 



/2/ M sinhe M V 
{ W ) 



Here we compare the imaginary parts of the phase shifts; Im(— iA) as: 



(5.2) 



Im(-«A string ) = 8Vp J. 



fc=l,3.5, 



k \ 4n 2 a'k 



__bf_ 2fc 

e 47ra' e 



n (i + /?) ■ 5 ft oath- (^). 



n=l 



(5.3) 



Im(-iA matrix ) = 8V P Y 



fc=l,3,5,- 



k 1 47r 2 a'fc 



" 4to' f M Binh e M 



n/* M 



Z=2 



A/ 



(5.4) 
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5.1 three special limits 

Now, if we assume a relation between the U(l) field strength F in string theory 
and the value of commutator F M among coordinates in IIB matrix model as 
follows: 



II (/r • i) H.// u * 



1=2 



1=2 



= 1 



considering three different regions of parameters which are b, v and f\ as: 

b 3> a' (long distance) 
v < 1 (low velocity) PJ, 
|/i| ^> f (large field strength) 



(5.5) 



(5.6) 



then, in each region, the both amplitudes approach the same value as follows: 



-r 3 



Vu 



1=2 



4na' 



(47r 2 a') 2 



1-ff 1-v 2 V b 2 



.(5.7) 



In addition , in the case of v <C 1 or |/i| 3> 1 the imaginary parts of the phase 
shifts; eq. (5.3) and eq. (5.4), are also identical. 



The main conclusion of this paper is that eq. (5.1) and eq. ( |5.2| ) are identical 
in three independent limit of parameters after identifying the field strengths in 
the models as eq (5.5). This condition is satisfied when we require 



det (r) + F)F 



M 



-1 



(5.8) 



We put fx pure imaginary which has the large absolute value, that is small real f£ 
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6 Conclusions and Discussions 



In this paper, we have considered two identical p-branes and compared their 
scattering amplitudes in the IIB string theory and the IIB matrix model. Here 
the p-brane has the lower brane charges in the various ways. That is, in the 
context of the string theory, there are arbitrary constant electromagnetic field 
strength F on the p-brane world volume. On the other hand, in the context 
of the matrix model, there are various values (F M ) of the commutators among 
brane coordinates. Based on the matrix model, we should calculate the am- 
plitude around the classical back-ground of the two body system of p-branes. 
However we have only the amplitude in the one loop level. The loop expansion 
is assumed to be the low energy expansion. Therefore the corresponding value 
based on the string theory is a cylinder diagram in the low energy limit, which 
is attached on the nontrivial background. It means the one loop calculation 
of open string with the corresponding boundary condition. After we have per- 
formed the low energy expansion that is the long distance expansion or small 
velocity expansion in each model, or small /* f expansion in IIB matrix model 
and large |/i| expansion in IIB string theory, the leading terms of the ampli- 
tudes are found to be the same precisely if we identify det (?/ + F)F M = —1. 
When all / ; M s are the same and fi = in the all world volume directions, our 
result reproduces that in Ref. ||. When / ; M <C 1, our result reproduces that in 
Ref [|[ [l2| . Note that we have performed these expansions of three parameters 
independently. In this paper, we have obtained a clear correspondence between 
the IIB string theory and its IIB matrix model in the near BPS configuration: 
D-brane scattering. Since we also find a way to study general boundary con- 
ditions we can calculate various configurations of two body systems of branes. 
For example, Dp-brane and Dg-brane case for p ^ q with different lower di- 
mensional brane charges and with arbitrary angles. We can also consider the 
relation between IIA string theory and matrix model for M-thcory fl3| . Our 
next task is to clarify the correspondence for objects which are far from BPS. If 
there is no agreement between string theory and matrix model for such cases, 
we may need to look for a better (matrix) model. 
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Appendix A 

In this appendix, we get commutation relations between modes' a m 's and x, 
namely, equations ( |2.25| )-( pi.27| ). The mode expansion for the coordinate shown 
in eq. (|2.2l|) is 



X(T,a) = x + MX{a + )+X(a-) 



X{p) 



-i{m-\-iE)p 



(A.l) 



To use eq. (A.l) not only in the region (0 < a < ir) but also in the region 
(— 7r < a < 7r), wc define extended variables in the following. 



f dX L (p) := MdX R {p) t-tt<p<t 
\dX R {p) := M- l dX L (p) t<p<t + tt ■ 

Then we get useful commutation relations as follows. 

[dX(p),dX(p')] = ma'5'(p - p')r) , r - tt < p,p' <r + tt 

A.l [a m ,a n ] 
Since we can write 

we get the commutation relation between a m 's as follows. 
1 2 rT+7T 



(A.2) 
(A.3) 



dpe i(m+iE)p dX(p) , 



(A.4) 



[a m , On] 



(2tt)2 a' J T _ n 
1 2 '• T+7r 



dpd// e 4(m+4jE)p aA(p), e iin+iE)p, dX(f/) 



, , . dpdp'e^+^ldXip^dXip')} (« 

— / T+ " dpd(/8'{p - p')e i(m+iE)p ^ f e^+^'V 
27T .L_^ V / 



i(n-\-iE)pf 



t—tt 

T + 7T 



dpdf/8(j> -(/)(?„ 



r\ e 



i(m+i_E)p\ I i(n+iE)p 



P—T+7T 



— I dp'S{p - p')e i( - m+iE ^r) f e l (" +iB ^' X ' 



(A.5) 



p—T — TT 



Since wc hnd that the last term in the eq. (A.5) is zero, the commutator is 



i 

"2^ 



dp (dp 



e i(m+iE)p | 



r — 7r 

2tt 



i(m+zi?)p 



i_B)p' 



(A.6) 
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The integrand in eq. (A. 6) is 



e i(m+n)p 



(A.7) 



because of M and M being SO(D — 1, 1) matrices. Therefore we get 

[otm, On] = S m+n (m + iE) 1] . (A. 8) 

A. 2 [x,a n ] 

We consider the following commutator 
[d T X(r,a),X(r,a')} = [dX(*-),x] 



M 
M 



dX{a + ),X(a + ) 



dX(a-),X(a- ) 



dX(a+),X(a- ) + dX(a-),X(a+ ) 



(A.9) 



Since 



[dX(p),X(p')] = ^e- i < m+iE )q am ,a„]| P «lp(e- i (» +i ^j 



= -mat (6(p- f/)-6(p))r) , 
wo rewrite the commutator as 



(A.10) 



[d T X{T,a),X{T,<r')} = [d T X( T ,a),x] 

- 2ma'5{(T - a')r) - ma' (M - M~ l ) 6 (a + a')?] 

+ nia' \s(a + ) (1 + M) + 5(a~) (l + M _1 )lrj . 



(ATI) 

Therefore 

[d T X(T,a),X(T,a')} = -2ma'S(a - a')r) , < a, a' < tt , (A.12) 



[8 t X(t, u), x] = —iria 

t 

[x, dX( P )] = ma'S(p) (1 + M) rj 



6{a+) (1 + M) + S(a~) (l + M~ l ) 



V 



(A.13) 



After integrating { e l ( m +' lE )p x eq. (A.13) } by p from r — n to r + 7r, we get 

(A.14) 
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A. 3 [x,x] 

We consider the following commutator 

[X(r,a),X( T ,a')] = [x,x] 

+ M[X(a + ),x] + [X(a-),x] 

x,X(a-')] + \x,X{a + ') 



M 
M 



A>+),A>+') 



M 



X(a-),X(a~ ) 



X(a + ),X(a-')] + \x(a-),X{*+') 



M 



(A.15) 



We find 

[X( P ),X(p')] 



and 



where 



TE I" d P I" dp'e-^ m +^P[a m ,a n ] ft 
1 m ,n Jo J ° V 

= ma' \e(p' - p) + e(p) - e(p') 
[x,X(p)] = nia' (1 + M) e(p)r) , 



— i(n J riE)p' 



'i - 



e(p):={0 P = 
"5 P <0 



(A.16) 
(A.17) 

(A.18) 



Therefore we get 

[X(t, <t),X(t, a')} = [x, x] - ina' (M - M' 1 ) e{o + a')rj = , < a, a' < n 



[x, x] = zf- (M - M- 1 ) j] 



(A.19) 



Appendix B 



In this appendix, we rewrite the commutation relations between the modes in 
a more convenient form and apply the formula to our case. Namely, we derive 
equations (|2.28|)-(|2.30|). We recall equations fl2.25|) -(|2.27D: 



[a m , a- 

[X : Ot-m 

\x,x] 



(m + iE) rj 
^f- (M - M~ l ) rj 



(B.l) 



2G 



We would like to find the conjugate pairs for the modes. Our strategy is to 



transform the each right-hand side of equations (B.l) to become the following 
form, respectively. 





m 







* 




(B.2) 



B.l T 

First, we would like to diagonalize the matrix E by transforming a m by a matrix 



E D 



:= T- l ET=: 



Then the commutation relations are 



E n 


Detail ^ 



\X,Or. 



(m + iE D ) j] T 5, 



m-\-n : 



i\l — (l + M) TrjT 



(B.3) 
(B.4) 
(B.5) 

(B.6) 
(B.7) 



B.2 S, N 



Second, we would like to write the each side of the equation (B.7) as the following 
form by transforming x 



* 

* * 



y := Sx . 

How to determine the matrix S and N is as the following. 



S := KJ 



K := 



-L 12 L 2 2 
1 



L:=JL, L:=(1 + M)T 



(B.8) 

(B.9) 
(B.10) 



if is a matrix which transforms L12 to 0, if L 22 exists. J is a matrix which 
exchanges the rows of a matrix (L) so that L 22 becomes existing 113. That is, 



17 The possibility is assumed here. If not so, we replace L22 by its part whose minor is not 
zero. In our case, we avoid pure Dirichlet directions: p + 2 • ■ ■ D — 1. 
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DetL 22 ^ 0. 
Then 

N := S{1 + M)T = KL 



Ln — i 12 i 22 1 L21 
L21 L 2 2 



N u 
N 2 i N 22 



Then the commutation relations are 



[y,a m ] = i\j—Nr) 



[y,y] 



B.3 C m --- step 1 



(NriT {ST) 1 - (ST) r]rN 



(B.ll) 

(B.12) 
(B.13) 



Third, we would like to write the each side of the equation (B.12) as the following 
form by transforming y. 





* 



Sm ■ 



We write y as follows: 



x := y 



(B.14) 



(B.15) 



Then 



\X. OL r , 



N u 

N 2 i N 22 



—m + iE\\ 
— m 



>]T 



Therefore we chose matrix C m as 



C m ^ := iN (m + iE D ) 



Co 



( N U E^ C12 

N 2 iE u 1 C22 



where cyi and C22 will be defined in the next subsection. 
Then we get 



\X,Ot. 



a' f 



2 10 JV : 



22 



(B.16) 



(B.17) 



(B.18) 
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B.4 C m --- step 2 



Finally, we would like to write the each side of the equation ( B.13 ) as the 
following form by choosing cyz and C22. 

^ o) ■ (B.19) 
We can find commutation relation for x's as follows: 

[x,x] = [y,y] 

Trig' ( Nu coth(7nEii)r7TiiiVH iVn coih^Eu^TnN^ 
2 \N 2 icoth.(TTEii)r] T iiNli N2icoth(irEu)riT n N^ 

+ M( _ci 2 iV| 2 

where we have used the formula 
x i 



N 22 c\ 2 c 22 Nl 2 -N 22 c\ 2 ) ' (B - 20) 
— 7r coth ttE\i . (B.21) 



' 12 



*-~L m + %En 

Therefore we choose Ci 2 and c 22 to satisfy following equations 

Cl2 iV* 2 = 7 riV 11 coth(7ri; 11 )?7 Tll iV* 1 +7r(^77 T (5r)*-(S*T)7 ?T ^* 
c 22 Nl 2 - N 22 ct 2 = 7riV 2 i cothfTr^O^n^x + 7T (iV^ST)* - (ST^tN*) 22 . 

(B.22) 

Then we get the following formula: 

TTIQ,' / Wii cothTrBiirjriiiV^ q 



[x,x] = y+NiinT 11 (ST)\^-{ST) linTll Nl 1 . (B.23) 

B.5 conjugate momentum for x 

We can define the conjugate momentum for the center of coordinate. 

[*>*d] = Vt(S ^ 22 ) • (R24) 

We assume that the determinant of the matrix i\r 2 2 is not zero. If not the case, 
we replace -/V 2 2 by its sub-matrix whose determinant is not zero. 
We define the momenta p l : 

p' := ^ (N^&oY , (B.25) 

where the number of the independent momenta is determined by the rank of 
N 22 . i is the index in the momentum space. Then 

[x\ P J ] = i8 13 .0 (B.26) 

18 If matrix S is not 77-othogonal matrix, we must find the momentum for x not x. 
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Now the trace of the zero mode in the phase space is defined as 

i 

where Li is a length of the area where X 1 covers M. 



(B.27) 



B.6 two identical parallel D-branes case 

We can find as follows after straight forward calculation 



m + ie 
m — ie 






ml 



o m +n , e := - hi 

7T \g-v 



iS lj , [x, a m7 ko] = [x,&q l ] = [x,x] = , 



1 + m|, •••,! + to'p+i p J , 



.9 := 



v / 2c7diagf-. 
yjl-f?(l-v 2 ) , i,j = 2---p + l, 



(B.28) 

where m-i, • • • , mp+2 are the 2x2 matrices which are defined in eq. (HH) and 



/i is a component of the field strength F in cq. ( p.3q ) 
0, 1, p + 1 components of matrices T, S, and iV are 



/ V2 



g+v 

g +v 



l 

V2 

1 

73g 



\ 



V2v 
g-v 

g-v 



9 ~) 





1 
1 
1 



2y/l-g 2 v 

g(g-v) 



(B.29) 



Note that [x, x] is zero as in Ref 



Appendix C 

In this appendix we find energy momentum tensor which satisfies the ordinary 



operator product expansion (OPE). Namely, we would like to explain eq. (3.21) 
and eq. ( 3.2 2| ) . For this purpose, we start with the calculation of OPE be 



tween naive energy momentum tensor of the world sheet. We normal order it 
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with respect to the vacuum defined in the section 2. We define a naive energy 
momentum tenser as follows. 



Tab := -47ra 



/ 1 ^^string 

7^ sh- b ' 



where String in (|C.l| ) is that of eq. (|3.l|) . 
With 



z := e 



we define the holomorphic energy momentum tensor: 
T(z) := ^(r o+T 01 ) 



1 



(ghost) 



(C.l) 



(C.2) 



^boson(^) -^fcrmion (^) ~t~ ^ghostt^) 

1 -m-n-2 T 

•= 2 Z : "m^" : 

+ ~ £ ^" s - 2 : finis + iE)d s : 
+ (ghost) . 

We extract the singular parts of the OPE in the following subsections. 



(C.3) 



C.l OPE for bosonic part 

In this subsection, we calculate OPE between the bosonic part of the naive 
energy momentum tensor as follows: 

T boson (z)T hoson {w) = J Tr (dZ (z)d Z (w)) 1 tit (dZ(z)dZ(w)) VT 

+ Tr (dZ(z)dZ(w)) t 7] T : dZ(z)dZ(w) : ry T 

+ regular( finite term with w — > z) , (C.4) 

where 



-^boson \%) 



dZ(z) T r]TdZ(z) 



-m—n—2 



(C.5) 



Z(z) := J-1*X{p = -i\*z) 
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E z~ n - lE °a n . 



First, we calculate the two point function for dZ: 

(dz( z y i dz(wy) . 

We calculate this separately as follows. 

(dZ{z)^Z{ W f) = -E(f) m ^(^-™) 



1 



(dZ(z) a dZ(w) b ) 



(z — w) 

1 



2 Vt 



zw 



^ OO 

-£(-) 

zw A — ' V z / 

m— 1 



where A, <5 = 2, • • • , L> - 1; a, 6 = 0, 1 and e a := (-) a e. 
In the case of (a, 6) = (0, 1), eq. flC.9|) is 



1 • / w \ m_ " oi 
> m - le) I - 1 r/ T 



m — 1 



ziu V z / 2 z — ' \ * 

5 

(?)' 



z 

yj \ ni—ie 



4 3 S 



r7 01 



In the case of (a, 6) = (1,0), eq. ( |C.9D is 

zw ^ 



If 

to + ie) I — 



m— 

- (-) d - E (- 

ziu V z J z * — ' V z 

m 

-)" 



Z 

w \ rn+ie 



Vt 



4»i 



T7 10 



That is, 



(dZ(zYdZ(w) v ) 



(HZ 
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(C.12) 

Therefore we get 

, , , . -j 2 / ie(l + ie) 

Tboson(Z)Tb oson {w) = 7 rj + 7 TJ Jboson(w) - 



ai 

(C.13) 



1 „ / ie(l + ie) . 
+ o w T hoson (w) — — — + regular , 

z — w \ 2w z 



where D = 10 is the space-time dimension. 
C.2 OPE for fermionic part 

We calculate OPE between the fermionic part of the energy momentum tensor: 

7f crmi on = -ttV : toWR ■ • (CM) 
2a' z 

We can calculate the two point function for ipR as follows. 

^ R (zr^ R (wY) = a'J2((z- {r+iE) dry (w- {s+iE) d s y) 

r,s 

r.s ^ a 



s<0 
r>0 



0<reZ-i^ ^ 

, r+iE 

Z 



0<reZ+i^- 



n=0 ' ' 



1 — — 2 

Z 

Therefore, after straight forward calculation, we get 

^"fermion (^)^fermion('^) 



V ■ (C15) 



D 
4 



(z — w) 
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+ 
+ 



[z — w) 2 
1 



-<l, 



z — w 
+ (regular) . 

Note that Tr E = is used in this calculation. 



^fermion ~l~ 
^fermion ~l~ 



aD -4TiE 2 



IQw 2 
ciD-ATy E 2 



IQw 2 



(C.16) 



C.3 shifted energy momentum tensor 

Summarizing previous two subsections, using Tr E 2 — 2e 2 , we define the follow- 
ing shifted energy momentum tensor: 

fii) := T(z)- ^\ + : €) + aD -^ 



2z 2 ' 16z 2 
1 ( aD ze N 
:- \ l(> 2 



(C.17) 



We can get the following OPE. 

f{z)f{w) = 



ID 



2T(w) + d w T{w) 



(z — w) 4 (z — w) 2 z — w 
+ regular + ghost contribution . 



(C.18) 



This form is the usual type and is invariant under the exchange of the z and w 
except for the regular part. 
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